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ABSTRACT 


This report presents a catalogue of theoretical geoidal 
variations for three types of structural features common to 
the earth's surface: seamounts , submarine ridges, and sub- 

marine trenches. These structures have been simulated by 
simple geometric shapes modeled in three-dimensions. A 
computer program calculated the potential and gravitational 
variations over the models. Profile plots of geoidal varia- 
tions and free-air gravity anomalies are presented over cross 
sections of the structures. The purpose of this catalogue is 
to provide a ready reference information set for comparison 
with satellite altimeter data for ocean areas. 
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INTRODUCTION 


The accurate delineation of large scale topography of the 
sea surface is equivalent to surface gravity measurements and 
hence can be useful for determining mass anomalies in the earth's 
crust. Satellite altimeter data offers two advantages over normal 
shipboard gravity observations. First, large scale features of 
several hundred to several thousand kilometer width are traversed 
quickly by a satellite pass, and hence problems of datum levels 
and instrument drift tend to be minimized. Secondly, orbiting 
satellites easily cross regions of the earth's surface that are 
remote and hazardous for surface travel. The southern oceans 
and the Arctic Ocean are particularly obvious regions in point. 

The radar altimeter information obtained during the SKYLAB mission 
instrument testing has already provided useful data on the marine 
geoid (McGoogan, et al, 1974; Leitao and McGoogan, in press; Hoge, 
1974? Bowin, 1974) . 

In order to assist in tho evaluation and interpretation of the 
SKYLAB radar altimeter data, as well as that from future missions, a 
catalogue of theoretical geoidal variations has been prepared for 
three types of structural features; seamounts, submarine ridges, 
and submarine trenches. The purpose of this catalogue is to pro- 
vide a ready reference information set for comparison with satellite 


altimeter data. 


METHOD 


For several years the gravity group at the Woods Hole Oceano- 
graphic Institution has been computing the gravitational attraction 
of models of crustal structure in their analyses of earth features. 
Recently a computer program (G3DC) was developed (Bowin et al, in 
press) to permit the rapid computation of the gravitational attrac- 
tion of three dimensional models which take into account the curva- 
ture of a planetary surface. This program was needed particularly 
for studies of the gravity field of the moon, because the moon has 
a much smaller radius than the earth and because gravity observa- 
tions are available at several heights above the lunar surface. 

Both these factors make it important to take into account the 
curvature of the moon’s surface. 

The potential of a mass varies inversely as the distance, 
whereas gravitational attraction varies inversely as the square of 
the distance. Accordingly, program logic was added to the three- 
dimensional gravity computation program mentioned above to compute 
the potential of each element of the model in addition to gravita- 
tional attraction. At each field point for which the computations 
are made, the potential Eind gravitational attraction are calcu- 
lated for each body of the model as well as for the complete model. 

If the densities used in the model calculations are those corresponding 
to mass anomalies from homogeneous spherical shells, then the potential 


values are those resulting from the anomalous mass and can be used 
directly in examining geoidal perturbations. If the full rock 
densities are used in the model, then the potential values obtained 
must be referenced to the potential value at some location where 
it is known or assumed. The results given in this catalogue were 
determined by the first method. 

Use of the computer program described above allows geoidal 
variations to be determined directly for any given model without 
the need to first compute a gravity field, and then from that 
field calculate geoidal variations . 

RESULTS 

The crustal structure models constructed for this catalogue 
are in complete isostatic equilibrium, since that is generally 
the condition of the earth's surface. For some of the models we 
present the potential and gravity anomalies separately for the 
topography and for its compensating mass at depth- In that way 
the effect of the topography can be seen as if the topography 
were a load (in the case of seamounts and ridges) or deficiency 
(in the case of trenches) on the earth's crust. Such results are 
labeled with a suffix T (for topography) or C (for compensation) - 
These results, of course, have larger magnitude variations of 
potential and gravity anomalies than do the models in isostatic 


equilibrium. 
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Three types of structures are presented: seamounts, simulated 

by vertical cones? ridges, simulated by horizontal triangular prisms 
and trenches, also simulated by horizontal triangular prisms, but 
with the apex pointing downward rather than upward. All of these 
models are constructed from polygonal laminae which stack one upon 
the other to comprise the bodies of the model. The technique of 
this type of model construction is described by Talwani and Ewing 
(1960) . If the laminae are thick, and the edges of the laminae 
are close to computation points, then the anomalies may show 
irregularities corresponding to the stair step edges of the three- 
dimensional model. 

Seamount models are labeled with the letter S, ridge models 
with the letter H, and trench models with the letter D. A summary 
tabulation in numeric order is given in Tables 4, 5, 6, 7, and 8. 

To help the reader more easily assimilate the results, three 
tables have been prepared to suggest some ways by which the models 
might be examined. The models shown in Table 1 help demonstrate 
the results for varying crustal density and height of the seamount. 
Ridge models having the same base width, crustal density, and 
height are indicated by the model numbers in parentheses. Thus 
the same table also aids an understanding of the difference in 
anomalies resulting from linear and circular features. 
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For the seamount models of Table 1, it is seen that the maxi- 
mum free-air gravity anomaly and the maximum geoid rise increases 
with increasing height of the seamount and increasing crustal 
density as would be expected. The same relations are true for 
the ridges. The magnitudes of these variations can be read from 
the model plots. A comparison of S-ll (water depth 4 km) with S-2 
(water depth 5 km) and S-12 with S-3 shows that decreasing the 
depth of the mass anomalies affects more the free-air anomaly than 
the potential (geoid) . 

For seamounts with a crustal density of 2.8 g/cm^ , a base 
width of 80 km, and water depth of 4 km, the following models 
show the effect of different seamount heights: S-10 (5 km height, 
island with 1 km relief) , S-ll (4 km height) , and S-12 (2 km height) . 
The magnitude of the gravity and geoid anomalies increases with 
increased height of the seamount. Note in S-10 that the anomalies 
plotted are for a line tangent to the shore of the island, not over 
the center of the feature. For seamounts with a constant height 
(5 km), water depth of 5 km, and a crustal density of 2.8 g/cm^, 
the following models show the effect of different base width: 3-15 
(40 km width), S-13 (60 km width), S-14 (100 km width), and S-16 
(120 km width) . As the width increases, the geoid rise increases 
in magnitude, but the free-air anomaly decreases. 

The effect of varying crustal density and base width for ridges 


having a constant height (4 km) in water 5 km deep can be studied 
by examining the models given in Table 2 . Free-air and geoid 
anomalies increase with increasing crustal density, but again, the 
free-air anomalies decrease whereas the geoid anomalies increase 
with increasing width of a ridge. Table 3 is an index to ridge 
models of constant crustal density (2.8 g/cm ) and constant water 
depth (5 km) but with varying ridge height and width. As the 
reader now has learned, the magnitude of both the anomalies in- 
creases with increasing height, but different relationships occur 
with increasing width. Although the maximum free-air anomaly value 
decreases with increasing width, the integral of the positive por- 
tion of the free-air anomaly curve increases with increasing width. 
Models with varying width but a different water depth (2 km) are 
also indexed in Table 3. 

The trench model numbers correspond directly with the radge 
model numbers except that heights for the ridges are depths for 
trenches. Thus Tables 1, 2, and 3 serve also as indexes to the 
trench models if the H prefix is changed to D. The trenches have 
a 25 km thick normal crust instead of 5 km as in the seamount and 
ridge models, in order to avoid having the compensation mass penetrate 
through the bottom of some of the trench models. The relations of 
free-air and geoid anomalies described previously apply also to the 
trench models but with the anomalies being negative instead of 
positive features. 
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The isostatically compensated trench model structure is not 
characteristic of the de^p-sea trenches associated with island arcs. 
These trenches are not in isostatic equilibrium, are sites of con- 
siderable mass deficiency, and hence are closer to the models with 
trench topography only (such as D-7T) . 

A comparison of the seamount and ridge models identified in 
Table 1 shows an interesting relation- For features of the same 
width and height, the ridge has a lower free-air gravity anomaly 
maximum but a higher geoidal rise than those of the seamount. 

This seemingly anomalous situation results from the greater gravi- 
tational effect of the compensation mass for the ridge than that for 
the seamount. It is only the vertical component of the gravita- 
tional attraction of anomalous masses in the model that contri- 
butes to the computed gravity anomaly values. Thus the greater 
effective mass deficiency of the compensation for ridge topography 
more than offsets the effect of its being at greater depth below 
sea level than the mass of the ridge topography. The free-air 
anomaly, therefore, is lower over the isostatically compensated 
ridge than over the isostatically compensated seamount of the 
same width even though the mass of the ridge is greater than that 


of the seamount. 
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Model 
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PARAMETERS FOR SEAMOUNT MODELS 
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PARAMETERS FOR TRENCH MODELS 
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